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Abstract 

Leptosphaeria maculans 'brassicae' is a damaging fungal pathogen of canola {Brassica napus), causing lesions on cotyledons 
and leaves, and cankers on the lower stem. A related species, L biglobosa 'canadensis', colonises cotyledons but causes few 
stem cankers. We describe the complement of genes encoding carbohydrate-active enzymes (CAZys) and peptidases of 
these fungi, as well as of four related plant pathogens. We also report dual-organism RNA-seq transcriptomes of these two 
Leptosphaeria species and B. napus during disease. During the first seven days of infection L. biglobosa 'canadensis', a 
necrotroph, expressed more cell wall degrading genes than L. maculans 'brassicae', a hemi-biotroph. L. maculans 'brassicae' 
expressed many genes in the Carbohydrate Binding Module class of CAZy, particularly CBM50 genes, with potential roles in 
the evasion of basal innate immunity in the host plant. At this time, three avirulence genes were amongst the top 20 most 
highly upregulated L. maculans 'brassicae' genes in planta. The two fungi had a similar number of peptidase genes, and 
trypsin was transcribed at high levels by both fungi early in infection. L. biglobosa 'canadensis' infection activated the 
jasmonic acid and salicylic acid defence pathways in B. napus, consistent with defence against necrotrophs. L. maculans 
'brassicae' triggered a high level of expression of isochorismate synthase 1 , a reporter for salicylic acid signalling. L. biglobosa 
'canadensis' infection triggered coordinated shutdown of photosynthesis genes, and a concomitant increase in transcription 
of cell wall remodelling genes of the host plant. Expression of particular classes of CAZy genes and the triggering of host 
defence and particular metabolic pathways are consistent with the necrotrophic lifestyle of L. biglobosa 'canadensis', and 
the hemibiotrophic life style of L. maculans 'brassicae'. 
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Introduction 

As more fungal genome sequences become available, it is 
apparent that their complement of genes and transcriptomes 
reflects fungal lifestyles. Lifestyles of plant pathogenic fungi are 
often classified into three broad categories: biotrophy, where the 
pathogen feeds from live host cells, necrotrophy, where the host 
cells are killed ahead of colonisation, and hemi-biotrophy, where 
the pathogen feeds from living cells before switching to a 
necrotrophic style of growth. These designations are imprecise 
and as the mechanisms of pathogenicity in a range of fungi are 
elucidated, lifestyle boundaries become more blurred [1]. 

The fungal genus Leptosphaeria belongs to the class Dothideo- 
mycetes, which includes a number of economically important 
plant pathogens that have a range of lifestyles on their hosts [2] . 
The Leptosphaeria species complex has two species L. maculans 



and L. biglobosa [3] and several sub-species or clades including L. 
maculans 'brassicae' and 'lepidii', and L. biglobosa 'canadensis', 
'brassicae', 'australensis' and 'occiaustralensis' [4]. The nomen- 
clature for these fungi is currendy under review [5] . L. maculans 
'brassicae', a hemibiotroph, causes blackleg, the most important 
disease of Brassica napus (canola) worldwide. Airborne sexual 
spores (ascospores) released from infested crop residues from the 
previous year's crop land on seedlings. Hyphae of germinated 
spores enter plant tissue via stomatal apertures and asymptomat- 
ically colonise the apoplastic spaces, between the plant cells. After 
eight to ten days, plant cells collapse and asexual sporulation 
begins within the necrotic leaf lesion. Hyphae then grow along the 
petiole and the stem, often resulting in a canker that girdles the 
stem causing lodging of the plant [6]. In contrast, L. biglobosa 
'canadensis' is not well-characterised. Although it causes cotyle- 
donary lesions, stem cankers are rarely produced [7]. 
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Leplosphaeria maculans 'brassicae' has numerous 'gene for 
gene' interactions with B. napus whereby an avirulence allele in 
the fungus renders it unable to attack cultivars with the 
corresponding resistance gene. This 'gene for gene' resistance is 
expressed in seedlings; five avirulence genes have been cloned so 
far - AvrLml , AvrLm4-7 , AvrLm6 AvrLml andAvrLmJl [8-12]. 
Only one resistance gene, LepR3, has been cloned from B. napus 
[13]. _ 

During infection fungi derive nutrition from the host plant, 
often by enzymatic degradation of proteins and carbohydrates. 
These latter enzymes are classed as Carbohydrate-Active enZymes 
(CAZys) and they often have well characterised domains. As well 
as providing nutrition, CAZy activity releases cell wall products 
that can act as DAMPs (Damage Associated Molecular Patterns) 
that activate the host immune system [14]. The biotrophic plant 
pathogens Blumeria graminis, Puccinia graminis, Melampsora 
laricus-populina, Hyaloperonospora arabidopsis and Ustilago 
maydis have fewer CAZy genes than necrotrophic plant pathogens 
do, perhaps because biotrophs do not need to digest plant cell 
walls for nutrition and they must evade the host immune system 
[15-17]. ^ 

Genomic sequences are now available for many Dothideomy- 
cetes and an extensive comparative analysis of 18 of them, 
including L. maculans 'brassicae' has been published [18]. These 
include the Brawica-infecting pathogen, Alternaria brassicicola, as 
well as three wheat-infecting pathogens, Stagonospora nodorum, 
Pyrenophora tritici-repentis, and Zymoseptoria tritici. The former 
three fungi, like many members of the order Pleosporales, have a 
necrotrophic lifestyle releasing toxins soon after invasion, whilst 
Zymoseptoria tritici, like most members of the order Capnodiales, 
is a hemibiotroph, with an extended period as a biotroph before 
causing necrosis [19-21]. 

The L. maculans 'brassicae' genome is compartmentalised into 
AT-rich -that are gene-poor comprising up to 35% of the genome, 
and gene- rich regions that have a high GC content [2], while the 
genome of L. biglobosa 'canadensis' is 30 Mb and lacks AT-rich 
regions (Grandaubert et al. manuscript submitted). A limited 
amount of oligo-array transcriptome data have been produced for 
L. maculans 'brassicae' during in vitro culture and infection of B. 
napus [1 1], but transcriptome data have not been reported for L. 
biglobosa or any closely related Brassica pathogen. 

Patterns of global gene expression can be generated by RNA- 
seq, a technique that enables analysis of dual transcriptomes; for 
instance, during a plant pathogen interaction. Furthermore RNA- 
seq can be exploited to analyse non-model organisms for which 
genomic resources are not well developed [22]. Few dual 
transcriptomes for plants and pathogenic fungi have been 
reported. Two recent reports of dual RNA-seq analysis of fungal 
diseases are of rice blast [23] and target leaf spot of sorghum [24]. 
Here we describe the genomes and transcriptomes of L. biglobosa 
'canadensis' and L. maculans 'brassicae' and canola [Brassica 
napus), during infection and in vitro. Our aim is to characterise 
genes that each pathogen uses to evade detection by the host, or to 
derive nutrition from the host, viz. the carbohydrate-active 
enzymes and peptidases. We also examine genes upregulated by 
the host during infection by each pathogen. 

Methods 

Fungal isolates and culture conditions 

L. maculans 'brassicae' isolate IBCN18 and L. biglobosa 
'canadensis' isolate 06J154, hereafter referred to as Lmb and 
Lbc, respectively, were subcultured on 10% Campbell's V8 juice 
agar at 22°C with a 12 h light/ 12 h dark light cycle. Conidia 



(5xl0 6 ) were added to 30 mL of liquid medium and incubated in 
still culture in a petri dish (15 cm diameter) at 22°C in the dark. 
Culture media were either 10% Campbell's V8 juice, or oilseed 
rape medium. The latter medium was prepared by homogenising 
leaves (200 g) of B. napus cv. Westar in a waring blender in a final 
volume of 1 L of water. The homogenate was centrifuged at 
2000 g for 20 min and the resulting supernatant was filter 
sterilised (0.22 um Millipore stericup filter). 

Plant growth and infection conditions 

Brassica napus cv. Westar was used for all infection assays; it 
has no known resistance genes. Seedlings were grown in a 
glasshouse maintained at 25°C under natural lighting. Wounded 
cotyledons were infected with conidia (10 uL of a 1 xlO' spores/ 
mL suspension) or water (mock inoculum), at 10 days post sowing 
as described previously [25]. 

Extraction of RNA and gene expression analysis 

For RNA-seq analysis, B. napus cotyledons were infected with 
Lmb, Lbc or water (control mock inocula), and at 7 and 14 days 
post inoculation (dpi) tissue around the inoculation site was 
harvested using a cork borer (0.5 cm diameter) and then placed 
into liquid nitrogen before freeze drying and subsequent grinding 
under liquid nitrogen. Tissue samples were prepared in biological 
triplicate. RNA was extracted using Trizol reagent from infected 
tissue and from mycelia of Lmb and Lbc from 7-day still cultures 
grown in oilseed rape medium. RNA was then DNAase-treated 
(Life Technologies) and cleaned up. 

The two biological replicates of each sample with the highest 
RNA integrity number values (>6) were sequenced with Illumina 
TruSeq version 3 chemistry on an Illumina HiSeq2000 sequencer 
at the Australian Genome Research Facility. In vitro derived RNA 
was sequenced with 1 00 bp paired-end reads in order to aid gene 
annotation, and in planta derived RNA was sequenced with 
100 bp single-end reads. A total of 15.5 Gbp sequence was 
generated from the in vitro libraries of the two fungi (7.75 Gbp per 
sample), and 72 Gbp sequence was generated from 12 in planta 
libraries (6 Gbp per sample) (Table SI in file S2). Reads were 
trimmed to a minimum phred quality score of 20 using Nesoni 
sequence software [26], orphaned members of pairs were retained, 
adaptor sequences were removed, and reads shorter than 20 bp 
were rejected. Trimmed reads were aligned to a reference genome 
sequence with Tophat vl.4. 1 splice-junction mapper [26]. 
Reference genomes were Lmb isolate v23.1.3 [2], Lbc isolate 
J154 (Grandaubert et al., submitted), and a Brassica exon array 
curated unigene set representing 135,201 gene models from B. 
napus, B. rapa and B. oleracea [27]. Aligned reads were quantified 
using Cufflinks vl.0.3 transcript assembly and quantification 
software and denoted as average expression levels (FPKM - 
fragments per kilobase of exon per million mapped reads) [28]. 
Cufflinks was run with bias correction to reduce variance due to 
sequence positional bias during Illumina sequencing. Gene 
expression FPKM values are listed in the relevant supplementary 
tables. All of the aligned RNA-seq reads have been deposited at 
the NCBI Sequence Read Archive (SRA), accessible under 
bioproject accession PRJNA230885 or sequence read archive 
SRP035525. SRA files may be read using the NCBI SRA toolkit 
(http:/ /www.ncbi.nlm.nih.gov/ sra). 

Quantitative RT-PCR experiments were carried out to 
determine levels of expression of genes containing CBM50 (LysM) 
domains. Cotyledon tissue from 32 B. napus cv. Westar seedlings 
was harvested at 3, 7 and 14 days after inoculation with Lmb 
isolate IBCN18. RNA (2 ug) was treated with DNase 1 (Life 
Technologies) for 1 h at 25°C, and reverse-transcribed using oligo- 
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Table 1. B. napus defence genes analysed [50]. 





Gene 


Full Name 


Defence Signalling Pathway 


B. nspus RRES unigene ID 


GenBank accession of CDS 


NCED3 


9-cis-epoxycarotenoid dioxygenase 3 


Abscisic acid 


rres040714.v1 


EV 137674 


ACS2a 


1 -amino - cyclopropanG - 1 -carboxylat6 
synthase 2 


Ethylene 


rresu/yoz/.v i 




CHI 


Chitinase 


Ethylene/salicylic acid//jasmonic acid 


rres036231.v1 


X61488 


HEL 


Hevein-like protein 


Ethylene/jasmonic acid 


rres036743.v1 


FG577475 


ICS1 


Isochorismate synthase 1 


Salicylic acid 


rres059693.v1 


EV225528 


PR-1 


Pathogenesis related protein 1 


Salicylic acid 


rres112514.v1 


BNU21849 


WRKY70 


WRKY transcription factor 70 


Salicylic acid 


rres038189.v1 


EV1 13862 


PDF1.2 


Plant Defensin 1.2 


Jasmonic acid 


rres071321.v1 


KC967203 



Defence signalling pathways indicate regulatory pathway to which the gene belongs. RRES unigenes are listed from the Brassica exon array [27], to which RNA-seq reads 
were aligned. Corresponding GenBank accession numbers are given for each gene. The response of each gene (with the exception of PDF1.2) to its corresponding 
signalling pathway was confirmed by Sasek et al [50]. 
doi:1 0.1 371 /journal.pone.01 03098.t001 



dT primer and Superscript III (Life Technologies) at 50°C for 1 h. 
Levels of gene expression were determined by qPCR using 
SensiMix (dT) SYBR Green PCR kit (Bioline) in a Corbett Rotor- 
Gene 3000 machine. Transcript levels of the gene of interest were 
normalized to that of Lmb actin as described previously (Gardiner 
et al., 2004). Primers are listed in Table S2 in file S2. 

Annotation of genes and domains 

Genes encoding CAZys were identified in six dothideomycetes 
(Lmb, Lbc, A. brassicicola, S. nodorum, P. tritici-repentis and Z. 
tritici) using www.cazy.org [29] and the dbCAN v3.0 HMM-based 
CAZy annotation server (http://csbl.bmb.uga.edu/dbCAN/) 
[30]. The major CAZy classes are Polysaccharide Lyases (PL), 
Glycosyl Transferases (GT), Glycosyl Hydrolases (GH), Carbohy- 
drate Esterases (CE), Carbohydrate Binding Modules (CBM) and 
Auxiliary Activities (AA). Secretion signal peptides were predicted 
using the SignalP 4.0 algorithm [31]. Pfam domains (Pfam A and 
B) were identified using profile hidden Markov models with 
HMMER3.0 [32] and Pfam_scan.pl software; the e value cut off 
was set to le-5 [33]. LysM-containing genes (with CBM50 
domain) were examined in more detail. Predicted gene models and 
aligned RNA-seq reads were viewed with the IGV browser [34], 
predicted intron splice sites were verified and translation start sites 
were checked for congruence with observed RNA-seq transcript 
boundaries. They were initially identified by comparison to the 
Pfam database [33], and then aligned to ECP6, the well 



characterised LysM-containing protein from Cladosporium fulvum 
[35]. Other Leptosphaeria proteins with LysM motifs were 
identified and characterised as described in supporting informa- 
tion (Figure SI in file SI). 

Peptidase domains were identified by BlastP [36] comparison of 
predicted protein sequences to the MEROPS 'peptidase database 
pepunit.lib dataset' of database of peptidase and inhibitor units 
(http://merops.sanger.ac.uk) [37] . 

Analysis of expression of fungal and B. napus genes 

The 100 most highly up regulated genes of Lmb and Lbc at 7 
and 14 dpi, compared to in vitro growth, were identified from the 
RNA-seq data (Table S 3-6 in file S2). Averages of quantile- 
normalised loglO-transformed FPKM values were calculated for 
each dataset for direct comparison of CAZy gene expression 
between Lmb and Lbc. Genes containing a CBM, GH, PL or AA 
domain and predicted to be secreted were analysed further. The 
sum of the expression values of each CAZy gene across the three 
treatments (in vitro, 7 and 14 dpi) was determined and the top 100 
genes were identified. Quantile-normalisation was applied so that 
expression values of genes of both fungi could be directly 
compared and loglO-transformed FPKM values were graphed. 
For each treatment, a heat map based on expression values was 
generated. 

B. napus defence genes, 9-cis-epoxycarotenoid dioxygenase 3 
(NCED), 1 -amino-cyclopropane- 1 -carboxylate synthase 2(ACS2), 



L. biglobosa 
J154 



L. maculans 
IBCN18 




3 dpi 5 dpi 7 dpi 10 dpi 14 dpi 17 dpi 



Figure 1. Symptoms on cotyledons of B. napus cv. Westar infected with L. maculans 'brassicae' or L. biglobosa 'canadensis'. B. napus cv. 
Westar cotyledons were wounded and inoculated with Lmb or Lbc spores and disease allowed to progress for 17 days post inoculation (dpi). 
Cotyledons were harvested and photographed at 3, 5, 7, 10, 14, and 17 dpi to track lesion development by the two pathogens. 
doi:1 0.1 371 /journal.pone.01 03098.g001 
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L. maculans L biglobosa Mock 
7 dpi 7 dpi 7 dpi 




Figure 2. Percentages of plant and fungal transcripts in B. napus cotyledons, uninoculated or at 7 and 14 days post-inoculation dpi 
with L maculans 'brassicae' (Lmb) or L biglobosa 'canadensis' (Lbc). B. napus cotyledons were infected with Lmb, Lbc or water (control mock 
inoculum). Total RNA was extracted from lesion tissue at 7 and 14 dpi for lllumina RNA-seq sequencing. The percentage of reads aligned to the 
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Brassica exon array unigene set (green) or the reference genomes for Lmb or Lbc (blue) is presented, as well as a photo of each representative 
infection. All of the aligned sequence reads were deposited at the NCBI Sequence Read Archive (SRA), accessible under bioproject accession 
SRP035525. 

doi:1 0.1 371 /journal.pone.01 03098.g002 



chitinase (CHI), hevein-like protein (HEL), isochorismate synthase 
1 (ICS1), Pathogenesis related protein 1(PR1), WRKY transcrip- 
tion factor 70, and plant defensin 1 (PDF 1-2) [38] were identified 
(Table 1). Their RNA-seq expression was analysed at 7 and 14 dpi 
in inoculated and uninoculated cotyledons. Additionally, expres- 
sion of Brassica genes involved in metabolic pathways was 
compared at 7 dpi by either Lmb or Lbc, and then analysed 
using MapMan, software that processes large gene expression 
datasets into metabolic pathways or other processes [39]. A 
Wilcoxon rank sum test with Bonferroni correction for multiple 
tests was used within MapMan to identify 20 functional categories 
of B. napus genes that were significantly regulated in response to 
infection by Lmb or Lbc. 

Results and Discussion 

Symptoms and lifestyles of L. biglobosa 'canadensis' and 
L. maculans 'brassicae' on cotyledons of B. napus 

Leptosphaeria maculans 'brassicae' (Lmb) and L. biglobosa 
'canadensis' (Lbc) exhibited different timing of symptom develop- 
ment on B. napus cotyledons (Figure 1). Lmb had a visually 
asymptomatic phase until after 7 days post inoculation (dpi), when 
lesions became visible. Lbc produced initial darkening of the 
cotyledon tissue at 3 dpi, and cell death and necrosis were 
apparent by 5 dpi and lesions increased in size until 1 7 dpi. At this 
time lesions caused by Lmb were of a similar size. Thus Lmb at 
7 dpi appeared to be growing biotrophically, but at 14 dpi was 
growing necrotrophically, whilst Lbc was necrotrophic from 5 days 
onwards, although as previously described, its growth is usually 
arrested before it can colonise the stem [7] . The rapid in planta 
necrosis caused by Lbc was similar to that previously described for 
L. biglobosa 'brassicae' [40]. 

General features of fungal and plant transcriptomes 

RNA-seq was used to define the transcriptomes of both 
pathogens and host during infection. A total of 87.5 Gbp of raw 
sequence was generated across all libraries (Table SI in file S2). In 
spite of the differences in disease symptoms at 7 dpi, the 
percentage of Lmb and Lbc reads aligned to the Lmb and Lbc 
reference genomes was similar (5 and 6% of the total), while 61 
and 44% of reads aligned to the reference genomes at 14 dpi 
(Figure 2). These data reflect that by 14 dpi the plant tissue is 
heavily colonised by both fungi. As expected, very few reads (less 
than 1%) from the mock-infected plant libraries aligned to the 
Leptosphaeria reference genomes. These libraries represent highly 
complex large datasets and Principal Component Analysis was 
carried out to compare the overall features of the transcriptomes, 
particularly gene identity and expression level. FPKM values were 
calculated for each gene in the three organisms. The duplicate sets 
of RNA-seq data were very similar, with the exception of the data 
for Lbc at 7 dpi. Gene expression of Lmb in the three conditions (7 
and 14 dpi, and in vitro growth) was clearly distinguishable one 
from another (Figure 3), and expression profiles of both fungi in 
vitro were more distinctive than in planta. Most of the variance 
(around 80%) was captured in PCI for Lmb (Figure 3A) and Lbc 
infections (Figure 3B), where the 7 and 14 dpi time points were 
distinguished from each other. The difference between in vitro and 
in planta samples was captured in PC2 in for both Lmb and Lbc 
plots. B. napus gene expression was markedly different at 7 dpi 



after infection by Lmb compared to Lbc, but similar at 14 dpi. 
The response of B. napus to Lmb at 7 dpi was similar to that of 
mock inoculation, implying that at this time the plant had not 
mounted a strong response to infection (Figure 3C). 

Since microarray expression (Nimblegen) data for another 
isolate of Lmb (v23.1.3) were available [2], the relative levels of 
expression of three avirulence genes (AvrLml, AvrLm4-7 and 
AvrLml 1), present in both Lmb isolates were compared to check 
for broad agreement between the two technologies (microarrays 
and RNA-seq). These avirulence genes were ranked in the top 100 
most highly expressed genes in planta for both microarray and 
RNA-seq, and both techniques reported a top 10 ranking for 
AvrLml and AvrLm4-7 (see below). 

Only eight of the 20 most highly upregulated in planta genes of 
Lmb had Pfam domains (Table 2). These included one gene with 
five CBM50 (LysM) domains (see later); other genes included three 
cytochrome P450 monooxygenases, and a transferase present in a 
gene cluster containing a polyketide synthase. Several hypothe- 
ticals were amongst the top 20 genes, as well as three avirulence 
genes, AvrLml, AvrLm4-7 and AvrLmJl [12]. In contrast, 19 of 
the 20 most highly upregulated in planta genes of Lbc had Pfam 
domains (Table 3), although five of these were conserved domains 
without functional annotation (Pfam-Bs, or Domain of Unknown 
Function (DUF)). Two cellulases, four other glycosyl hydrolases, as 
well as three peptidases were present. 

At 14 dpi, 14 of the top 20 in planta-expressed genes of Lmb 
were hypothetical genes lacking any functional annotation 
(Table 4). Few were also in the top 20 at 7 dpi, with the exception 
of AvrLml and LemaPl 14790.1, which has no Pfam domains and 
was the most highly expressed gene at both time points. Even 
though there was a high degree of necrosis on cotyledons at 
1 4 dpi, only two hydrolytic enzymes (trypsin and a glycohydrolase 
7) were included in the top 20 most highly expressed Lmb genes. 
Two of the top 20 genes of Lbc were CAZys; two peptidases and 
three dehydrogenases, genes that may have degrading roles were 
also present (Table 5). 

Small secreted proteins (SSPs) were well represented in the most 
highly in planta upregulated genes of Lmb. These included 
hypothetical, avirulence and CAZy genes. Ohm et al (2012) 
compared the small secreted proteins of 1 8 Dothideomycete fungi 
including Lmb isolate 23.1.3, and found that 21.3% of all Lmb 
SSPs identified were singletons with no homologue in any of the 
other genomes [18]. A recent study found 30 of the 100 most 
highly expressed SSPs in planta at 7 dpi were unique to Lmb [52]. 

In summary, amongst the top 20 most highly expressed genes in 
planta, fewer encoding cell wall degrading enzymes were 
expressed by Lmb than by Lbc at 7 and 14 dpi. Cell wall 
degrading enzymes not only provide nutrition for the pathogen by 
hydrolysing carbohydrates, but facilitate fungal progression 
through the plant apoplast (intercellular spaces), during biotrophic 
stages of infection. The repertoire and expression of fungal CAZy 
genes presumably reflects the cognate carbohydrate present in the 
host plant and it is notable that enzymes such as those degrading 
cellulose are very highly expressed. The cell wall of leaves of 
Arabidopsis thaliana, which like B. napus is a crucifer, contain 
polysaccharides including rhamnogalacturonan I and II (as 
pectin), xyloglucan, glucuronoarabinoxylan, and cellulose (14%). 
Presumably the cotyledon has a similar polysaccharide profile. An 
additional 14% of the wall is composed of protein [41]. More than 



PLOS ONE | www.plosone.org 



5 



July 2014 | Volume 9 | Issue 7 | e103098 



O 

O 
0- 



B 



CM 

o 

Q_ 



CM 

o 

CL 



1.0 



0.5 - 



0.0 



-0.5 



-1.0 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 
PC1 (83 %) 




0 0.1 0.2 0.3 0.4 0.5 0.6 
PC1 (86 %) 




0.2 0.3 0.4 0.5 0.6 



PC1 (86 %) 



PLOS ONE | www.plosone.org 



The Leptosphaeria Transcriptome 



Figure 3. Principal Component Analyses of duplicate sets of 
RNA-seq data for L maculans 'brassicae' (Lmb), L. biglobosa 
'canadensis' (Lbc) and B. napus. Scores plots for principal 
component analysis of expression values (FPKM) of genes of Lmb 
IBCN18 (A), Lbc J154 (B), or B. napus (C). Gene expression values during 
infections of cotyledons are plotted in orange (7 dpi) and green 
(14 dpi). Gene expression values during growth in vitro are plotted in 
purple; those in mock-infected cotyledons are plotted in magenta. The 
percentage of the total variance explained is listed on each axis label. 
doi:10.1 371/journal.pone.01 03098.g003 

a third of the carbohydrate is soluble in phosphate buffered saline, 
including half of the total pectin, suggesting it would be readily 
available to a pathogen growing in the apoplast [41]. 

Distribution of CAZy domains and their expression 

In view of the dominance of CAZys in the 20 most highly 
expressed genes of Lbc after seven days in planta, CAZy domains 
were sought in genome sequences of Lbc, Lmb, and also in four 
other Dothideomycetes for comparison. The repertoire of domains 
derived from the dbCAN database and its HMM-based sequence 
similarity search was consistent with CAZy classifications carried 
out previously on Dothideomycetes (Table S7 in file S2) [18,42], 
and the www.cazy.org annotation for Lmb isolate v23.1.3. In 
general, the numbers of carbohydrate-binding module (CBM), 
glycosyl hydrolase (GH), glycosyl transferase (GT), and polysac- 
charide lyase (PL) domains were in agreement with previous 
studies (Table 6). CAZys are divided into sub classes, generally 
based on substrate specificity or ligand; some of these are 
described later. Two CAZy classes in Lmb were identified more 
frequently using our analyses; CE domains (109 vs 34) and AA (78 
vs 27). The difference in numbers of CE domains was probably 
because 27 CE1 and 42 CE10 domains had been previously 
excluded from the CAZy expert curated dataset on the basis that 
they were likely to act upon non-carbohydrate substrates [43] . AA 
domains are a newly formed CAZy classification for ligninolytic or 
lytic polysaccharide monooxygenases. Our automated approach 
may have included closely related monooxygenases that do not 
degrade lignin or polysaccharides. 

Most of the CAZy containing-genes had only one CAZy 
domain, but several had multiple CBM domains, particularly of 
the subclass CBM50 (see later and Figure SI in file SI). Five of the 
six dothideomycetes had between 580 and 700 CAZy domains, 
but Z. tritici only had 489, which is in general agreement with 
previous findings [44] . Overall the two Leptosphaeria species had 
similar complements of CAZy genes. Alternaria brassicicola had 
the next most similar profile. Pyrenophora tritici-repentis had a 
profile more similar to Stagonospora nodorum than to Lepto- 
sphaeria, consistent with the host of the two former fungi being a 
monocotyledonous cereal, rather than a dicotyledonous oilseed 
plant. Dicot cell walls generally contain a higher proportion of 
pectin, compared to glucuronoarabinoxylan, which is predomi- 
nant in a typical monocot cell wall and believed to partially 
substitute for low levels of pectin in monocot cell walls [15]. 
Zymoseptoria tritici had a distinct profile with many fewer CAZy 
genes. This is congruent with previously reported characteristics of 
the members of order Capnodiales, versus the Pleosporales, to 
which the other five fungi belong. L maculans 'brassicae' had 
fewer AA domains than Lbc, S. nodorum and P. tritici-repentis 
did. The number of PL domains ranged from 19 to 24 in Lmb, 
Lbc and A. brassicicola. In contrast, the three wheat pathogens, S. 
nodorum, P. tritici repentis and Z. tritici had many fewer - 
between 4 and 10. 

Although the complement of CAZys in the two Leptosphaeria 
species was generally similar, there were some interesting 
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Table 6. Distribution of CAZy domains in predicted proteins of six Dothideomycetes. 





Number of domains 


CAZy domain 


AA 


CBM 


CE 


GH 


GT 


PL 


Total 


Fungus (lifestyle) 


Leptosphaeria biglobosa (necrotroph) 


98 


85 


117 


234 


106 


23 


663 


Leptosphaeria maculans (hemibiotroph) 


78 


63 


109 


217 


100 


19 


586 


Alternaria brassicicola (necrotroph) 


92 


66 


120 


233 


92 


24 


627 


Stagonospora nodorum (necrotroph) 


122 


66 


143 


267 


96 


10 


704 


Pyrenophora tritici-repentis (necrotroph) 


114 


56 


124 


246 


105 


10 


655 


Zymoseptoria tritici (hemibiotroph) 


55 


26 


97 


199 


108 


4 


489 



The classes of CAZy domains are Auxiliary Activities (AA) Carbohydrate Binding Modules (CBM), Carbohydrate Esterases (CE), Glycosyl Hydrolases (GH), Glycosyl 
Tranferases (GT) Polysaccharide lyases (PL). Some genes contain multiple CAZy domains. 
doi:1 0.1 371 /joumal.pone.01 03098.t006 



differences. The major difference was in the CBM class with > 
30% more domains present in Lbc than in Lmb. L maculans 
'brassicae' had fewer CBM18-containing-genes than Lbc, but a 
similar number of CBM 18 domains. CBM 18 domains in Lmb 
were in a more compact gene structure (28 domains in 1 1 genes); 
indeed this domain was always present in multiple copies within a 
gene, five genes contained homopolymers of only CBM18. L. 
biglobosa 'canadensis' had 32 domains in 21 genes, eight of which 
had a single CBM 18 domain. CBM50 (LysM), CE4, and GH18 
CAZy domains occurred frequently in genes that had CBM 18 
domains. (Table S8 in file S2). In eukaryotes this domain often has 
a chitm-binding function and modules are often adjacent to 
chitinase catalytic domains, but this domain is also present in non- 
catalytic proteins either singly or as multiple repeats. CBM 18 
domains may be involved in targeting the degradative domains to 
particular carbohydrates, or may enhance catalysis by ensuring 
retention of the substrate between catalytic cycles. 

About half of the genes with CAZy domains in the two 
Leptosphaeria species did not have a signal peptide and thus were 
not predicted as being secreted. This may be due to incorrect 
annotation at the 5' end of genes, genes incorrectly merged during 
annotation, or a false negative result from SignalP. Pectate lyases 
were the most frequently secreted CAZy (100% for Lbc, 89% for 
Lmb), while GTs were the least frequently secreted (1 1 % for Lbc, 
12% for Lmb). This is consistent with the role of GTs in 
synthesising carbohydrates for the fungal cell wall, and their 
membrane or intracellular location. 



We then examined the RNA-seq data to determine expression 
levels of CAZys of Lbc and Lmb during infection of B. napus 
cotyledons (Table 7, Table S7 in file S2). The expression profiles 
of the top 100 most highly expressed (summed FPKM values at 
7 dpi, 14 dpi, and in vitro) secreted CAZy genes are presented in 
Figure 4. At 7 dpi, Lmb had a high level of expression of genes 
with CBM domains, particularly CBM50/LysM (for instance, 
genes A and C in Figure 4); this level decreased by more than 50% 
by 14 dpi and was even lower during growth in vitro. 

Genes with CBM domains were also among the most highly 
expressed CAZy genes of Lbc in planta, but the overall level of 
expression was much lower than that of Lmb genes, at the same 
time point. To reveal the most in planta specific CAZy classes, 
expression ratios between in planta and in vitro growth were 
calculated for each CAZy class (Table S9 in file S2). For Lmb, the 
CBM50 (chitin-binding) domain was highly upregulated in planta 
at 7 dpi, whereas at the same time point Lbc upregulated CBM6 
(cellulose binding) domains more highly (Table S9 in file S2). The 
in planta upregulated RNA-seq expression profile of the CBM50- 
containing genes of Lmb was validated by quantitative RT-PCR 
(Figure S2 in file SI). These experiments showed peak expression 
of each LysM gene at 7 dpi in planta and much lower expression 
levels in vitro. At 3 dpi, LysM gene expression was higher than in 
vitro levels, but lower than levels at 7 dpi. For the CBM 18- 
containing genes, expression was varied, with >3 orders of 
magnitude between the highest and lowest expressed genes. L. 
biglobosa 'canadensis' had generally higher expression compared 



Table 7. Expression of CAZy genes in L, maculans 'brassicae' and L biglobosa 'canadensis.' 







L. maculans 






L biglobosa 






CAZy class 


7dpi in planta (FPKM) 


14dpi in planta (FPKM) 


in vitro (FPKM) 


7dpi in planta (FPKM) 


14dpi in planta (FPKM) 


in vitro (FPKM) 


AA 


86.5 


109.9 


117.3 


104.5 


89.1 


159.7 


CBM 


508.7 


203.6 


158.5 


163.7 


174.3 


341.4 


CE 


78.1 


69.5 


63.6 


78.7 


55.8 


99.9 


GH 


110.2 


122.3 


88.1 


122.7 


81.3 


125.5 


GT 


54.8 


60.9 


57.1 


46.3 


50.0 


112.3 


PL 


56.3 


23.9 


6.9 


121.3 


37.7 


29.4 



RNA-seq derived gene expression values were mapped onto the identified CAZy domains for Lmb and Lbc. Average expression levels (FPKM) were then calculated for 

each major class of CAZy. 

doi:1 0.1 371 /joumal.pone.01 03098.t007 
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in vitro 
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in vitro 
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in planta 
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and in vitro 



PLOS ONE | www.plosone.org 



12 



July 2014 | Volume 9 | Issue 7 | e103098 



The Leptosphaeria Transcriptome 



Figure 4. Expression profiles of secreted CAZys of L. maculans 'brassicae' (Lmb), L. biglobosa 'canadensis' (Lbc) in planta and in vitro. 
The top 100 genes expressed across the three treatments (7, 14 dpi and in vitro) and predicted to be secreted and to contain a CAZy domain (CBM, 
GH, PL, or AA) were selected. Quantile-normalisation was applied and log1 0-transformed FPKM values were graphed. The intensity of blue shading is 
proportional to the expression level. The gene order is based on a dendrogram created from a Euclidean similarity matrix with average group 
distance. Letters with a triangle point to genes that were amongst the top 20 most highly expressed genes in vitro or in planta (Tables 2, 3, 4, or 5). A/ 
Lema_P1 02640.1 (Lm2LysM); B/Lema_P01 3700.1 (Glycoside hydrolase family 7); C/Lema_P070100.1 (Lm5LysM); D/LbJ154_P005286 (pectate lyase); 
E/LbJ154_P009204 (cellulase); F/LbJ154_P001246 (Glycosyl hydrolase family 43); G/LbJ154_P001652 (cellulase); H/LbJ154_P009247 (Glycosyl 
hydrolases family 12); l/LbJ154_P004093 Glycosyl hydrolases family 39). Categories of gene expression (High expression in vitro, High expression in 
planta, Medium expression in planta and in vitro, High expression in planta at 14 dpi, High expression in planta at 7dpi, High expression in planta and 
in vitro) were manually assigned and indicated by coloured polygons linking the genes in each category across Lmb and Lbc; the number of genes in 
each category is indicated on the vertical sides of the polygon. 
doi:10.1371/journal.pone.0103098.g004 



to that of Lmb (Table S8 in file S2). For Lbc, the average 
expression value of all CBM genes was highest during in vitro 
growth, which was due to high expression of several genes that had 
both CBMs and hydrolytic CAZY domains (Table 7). Such genes 
included Lb_j 154_P004089, which has three CBM 18 and one 
CE4 domains (Table S8 in file S2). 

Polysaccharide lyases (PL) were expressed highly at 7 dpi, and 
expression decreased by 14 dpi in both species. At 7 dpi, the 
expression levels of PLs in Lbc were twice those of Lmb. 
Expression levels then dropped by 75% at 14 dpi and were low 
during growth in vitro. L. biglobosa 'canadensis' had higher levels 
of expression of degrading CAZys such as AAs, GHs and PLs 
(genes D,E,F,G,H and I; Figure 4) at 7 dpi, than at 14 dpi, when 
lesion growth had slowed. The PL3 class (pectate lyase) was highly 
upregulated by Lbc at 7 dpi (Table S9 in file S2), and was also the 
most upregulated PL class in Lmb. This may be a reflection of the 
abundance of pectin in the cotyledon. 

Glycosyl transferase (GT) -containing genes were expressed at 
similar levels in all conditions in both Leptosphaeria species, except 
for a two-fold increase in Lbc grown in vitro compared to in 
planta (Table 7). The class GT21, which encodes biosynthetic 
enzymes that glycosylate lipids, was the most in planta upregulated 
class of glycosyltransferases. Of the carbohydrate esterase 
enzymes, CE8 (pectin methylesterase) and CE12 (pectin acetyles- 
terase) were the most in planta upregulated classes for Lbc and 
Lmb, respectively (Table S9 in file S2). 

At 1 4 dpi, of the six CAZy classes, expression of only PLs was 
higher in Lbc than in Lmb. By 14 dpi, expression of CBM and PL- 
containing genes by Lmb had decreased, but expression of genes 
with AA and GH domains had increased. In vitro growth was 
characterised by low CBM expression and very low PL expression 
(Table 7). The CAZy expression profile for Lmb in planta may 
reflect early avoidance of chitin- triggered immunity via high 
expression of CBMs such as LysM genes, whilst the fungus feeds 
from soluble pectin via PL genes. At 14 dpi elevated expression of 



the oxidative AA class and the hydrolytic GH class may facilitate 
degradation of lignin and cellulose in dead cells. 

In many fungi, expression of cell wall degrading enzymes 
CAZys is regulated by carbon catabolite repression, a global 
mechanism that ensures readily assimilated carbon sources such as 
glucose are preferentially used. A Lmb homologue of the 
Saccharomyces cerevisiae sucrose non-fermenting protein kinase 1 
C SNF1 ) carbon catabolite regulator has been recently shown to 
regulate expression of several CAZy genes encoding pectate lyases, 
beta-l,3-glucanase, and a glucosidase [45]. These CAZy genes 
were upregulated four days after inoculation of canola cotyledons 
with a wild type strain of Lmb; whereas an LmSNFl knockout 
strain had significandy reduced expression of these CAZys during 
growth on pectin in vitro. We recorded similar upregulation of the 
CAZy genes encoding two pectate lyases and the carbohydrate 
esterase in planta at 7 dpi, but not significant upregulation of the 
chitin deacetylase. This latter difference could be due to differing 
time points at which tissue was analysed and the more robust 
normalisation used in an RJNA-seq analysis compared to a qPCR 
method. 

In both Lmb and Lbc, the top 5% most highly expressed CAZy 
genes accounted for approximately 50% of the sum total CAZy 
expression (FPKM) at each growth stage (7 and 1 4 dpi and in vitro 
growth) that was analysed (data not shown). This suggests that the 
vast majority of CAZys are expressed at low levels or only turned 
on at specific times in the fungal life cycle. This phenomenon has 
been reported previously for Z. tritici, which differentially 
expresses CAZy genes, such as cutinases (CE5 subclass) according 
to biotrophic, necrotrophic or saprotrophic stages of growth [44]. 
The hemi-biotroph Collelotrichum higginsianum upregulates 
expression of effectors and secondary metabolite biosynthetic 
enzymes both before penetration and during biotrophic growth on 
Arabidopsis. C. higginsaneum then upregulates hydrolases and 
transporters at a later stage, each wave delivered according to the 
stage of pathogenic transition [46]. This again matches our 



Table 8. Distribution of peptidase domains in six Dothideomycete genomes. 



Peptidase type 


A 


C 


G 1 


M 


S 


T U Total 


Species (Lifestyle) 


Leptosphaeria biglobosa (necrotroph) 


18 


73 


1 8 


144 


145 


23 1 413 


Leptosphaeria maculans (hemibiotroph) 


16 


78 


1 10 


140 


140 


25 1 411 


Alternaria brassicicola (necrotroph) 


17 


77 


1 6 


146 


163 


23 1 434 


Stagonospora nodorum (necrotroph) 


23 


80 


0 7 


151 


192 


23 1 477 


Pyrenophora tritici-repentis (necrotroph) 


22 


78 


1 9 


155 


163 


22 1 451 


Zymoseptoria tritici (hemibiotroph) 


33 


67 


4 8 


130 


191 


23 1 457 


Peptidase domains are classified according to the 


catalytic type, or 


nhibitor activity. Peptidase types are Aspartic (A), Cysteine (C), Glutamic (G), Inhibitor, Mettallo-(M), 



Serine (S), Threonine (T), Unknown (U). 
doi:1 0.1 371 /journal.pone.01 03098.W08 
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Figure 5. Expression of key defence genes of B. napus at 7 or 14 dpi with either L. maculans 'brassicae' (Lmb) or L. biglobosa 
'canadensis' (Lbc). RNA-seq data for eight key Brassica defence genes were determined at 7 and 14 dpi. Average expression (FPKM) after infection 
by L maculans (blue), L biglobosa (red), and mock inoculum (green) is plotted. Error bars indicate the Cufflinks 95% confidence interval for each FPKM 
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value. The genes assayed were 9-cis-epoxycarotenoid dioxygenase 3 (NCED3), 1 -amino-cyclopropane-1 -carboxylate synthase 2 (ACS2), chitinase (CHI), 
hevein-like protein (HEL), isochorismate synthase 1 (ICS1), Pathogenesis related protein 1 (PR-1), WRKY transcription factor 70 (WRKY70), and plant 
defensin 1 (PDF1.2). Each gene was also classified according to hormone(s) abscisic acid (ABA), ethylene (ET), jasmonic acid (JA) and salicylic acid (SA) 
that induced higher expression. 
doi:1 0.1 371 /journal.pone.01 03098.g005 



observation that L. biglobosa lacks or only has a short biotrophic 
stage and expresses hydrolases earlier than Lmb does. CAZy genes 
of the two Leptosphaeria species that have low expression during 
development of cotyledonary lesions might be expressed more 
highly in another growth situation, such as earlier in penetration of 
the tissue, or during colonisation of petiole, stem, or during 
saprophytic growth on woody stubble of B. napus. 

Peptidase distribution and expression 

Since the two Leptosphaeria species grow between the plant cells 
in the apoplast, protein from the plant cell wall is a potential 
source of nitrogen for nutrition. A trypsin-like peptidase was highly 
expressed during necrotic stages of infection for both Lmb (14 dpi) 
and Lbc (7 dpi) (Table 2, Table 3), therefore we characterised the 
peptidase content of both Leptosphaeria genomes and compared 
them to four other Dothideomycete plant pathogens. 

The MEROPS peptidase classifications categorises peptidases 
by their catalytic nucleophile [37]. For example, aspartic (A), 
cysteine (C), serine (S), threonine (T) peptidases are named for the 
corresponding catalytic residue. Exceptions are for metallo (M) 



peptidases, which use a coordinated metal ion, and the classifi- 
cations for unknown activities (U) and peptidase inhibitors (I). 
Peptidase composition and number was very similar for Lmb, Lbc 
and A. brassicicola (Table 8). Metallo- and serine-peptidase genes 
comprised almost 70% of the total peptidases. The next most 
abundant classes were the cysteine, threonine and aspartic 
peptidases, in that order. L. biglobosa 'canadensis' had seven 
C56 cysteine peptidases (7) whilst Lmb only had three. C56 
peptidases act on peptides <20 amino acids, presumably requiring 
prior digestion of the substrate by another peptidase. P. tritici- 
repentis and S. nodorum had more peptidases in the metallo and 
serine classes than Lmb or Lbc. The hemi-biotroph Z. tritici had a 
larger complement of serine and aspartic peptidases than the two 
Leptosphaeria species did. 

In general, the expression levels of peptidase at 7 and 1 4 dpi in 
planta were similar in Lmb and Lbc (Pearson correlation 
coefficient of 0.76 for Lmb IBCN18 and 0.85 for LbcJ154), with 
a few notable exceptions. A trypsin-like serine peptidase (MER- 
OPS S01 A) was the most highly in planta up-regulated peptidase 
(500-fold higher at 7 dpi in planta than in vitro) in both Lbc 




-3 0 +3 

L. biglobosa 1 I L. maculans 

Brassica napus gene expression 
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Figure 6. Transcription of B. napus genes involved in metabolic processes including photosynthesis seven days after inoculation 
with L. maculans 'brassicae' (Lmb) or L biglobosa 'canadensis' (Lbc). RNA-seq gene expression values for a 8. napus unigene set [27] were 
used to calculate a ratio of expression values (log2) for 8. napus genes after infection by Lmb or Lbc. Ratios were plotted on major metabolic 
pathways with Mapman software [39]. A yellow square indicates a 8. napus gene that is expressed more highly during Lmb infection, while a blue 
square indicates a 8. napus gene with higher expression during Lbc infection. An expression ratio close to zero is shown with a white square and 
indicates equivalent expression during infection by either pathogen. Only genes with expression values greater than 10 FPKM were included. 
Abbreviations: LDH, lactate dehydrogenase; ADH, Alcohol dehydrogenases, TCA, tricarboxylic acid cycle; raff, raffinose; Treh, trehalose; PSI, 
photosystem one; PSII, photosystem two; ABA, abscisic acid; ET, Ethylene; SA, salicylic acid; JA, jasmonic acid. 
doi:1 0.1 371 /journal.pone.01 03098.g006 
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Table 9. Top 20 functional categories that are significantly regulated in B. napus in response to infection by L maculans 'brassicae' 
compared to infection by L. biglobosa 'canadensis.' 



Mapman Function 
Category 


Mapman Function Description 


Number of 
genes 


p-value 


1 


Photosynthesis 


313 


0.00E+00 


1.1 


Photosynthesis, light reactions 


210 


0.00E+00 


1.1.1 


Photosynthesis, light reactions, photosystem II 


105 


0.00E+00 


26 


Miscellaneous 


365 


4.42E-25 


10 


Cell wall 


116 


4.31 E-1 2 


1.1.2 


Photosynthesis lightreaction photosystem 1 


44 


5.06E-1 1 


1.1.1.2 


Photosynthesis lightreaction. photosystem II, PS1I polypeptide subunits 


63 


5.06E-1 1 


1.3 


Photosynthesis calvin cycle 


76 


2.19E-10 


1.1.1.1 


PS lightreactions photosystem II, LHC-II 


42 


3.08E-10 


1.1.2.2 


PS lightreactions photosystem 1, PSI polypeptide subunits 


39 


3.56E-10 


20.1 


Stress, biotic 


123 


1 .98E-09 


10.6 


Cell wall, degradation 


28 


3.88E-09 


26.12 


Misc peroxidases 


26 


1.07E-08 


26.8 


Misc nitrilases, nitrile lyases, berberine bridge enzymes, reticuline oxidases, troponine reductases 


42 


3.72E-08 


26.10 


Misc cytochrome P450 


46 


1.93E-07 


26.3 


Misc gluco-, galacto- and mannosidases 


42 


1.35E-06 


10.6.3 


Cell wall degradation, pectate lyases and polygalacturonases 


17 


1 .74E-05 


35 


Not assigned 


1969 


6.69E-05 


35.2 


Not assigned, unknown 


1969 


6.69E-05 


17 


Hormone metabolism 


107 


8.21 E-05 


2.2.1 


Major CHO metabolism, degradation sucrose 


22 


6.47 E-04 


20.1.7.6 


Stress biotic, PR-proteins, proteinase inhibitors 


12 


6.63 E-04 


1.2 


Photosynthesis, photorespiration 


27 


7.32E-04 


17.5 


Hormone metabolism.ethylene 


30 


7.77E-04 


34 


Transport 


208 


8.58E-04 


33 


Development 


107 


1.62E-03 


26.9 


Misc glutathione-S-transferases 


32 


1.92E-03 


10.8 


Cell wall pectin esterases 


18 


1.92E-03 


17.5.1 


Hormone metabolism, ethylene synthesis-degradation 


21 


1.92E-03 


20.1.7.6.1 


Stress biotic, PR-proteins, proteinase inhibitors, trypsin inhibitor 


11 


1.92E-03 


1.1.3 


Photosynthesis, light reactions, cytochrome b6/f 


8 


2.25E-03 


10.8.1 


Cell wall, pectin esterases, PME 


13 


2.35E-03 


31.4 


Cell vesicle transport 


32 


2.41 E-03 


1.3.6 


Photosynthesis, calvin cycle,aldolase 


13 


2.49E-03 


1.1.4 


Photosynthesis, lightreaction, ATP synthase 


15 


2.55E-03 


29.5.3 


Protein degradation, cysteine protease 


18 


3.1 1 E-03 


29.5.9 


Protein degradation, AAA type 


15 


3.41 E-03 


2.2.1.3 


Major CHO metabolism, degradation, sucrose invertases 


13 


4.18E-03 


13.2.3 


Amino acid metabolism, degradation, aspartate family 


11 


5.62E-03 


1.1.5 


Photosynthesis, light reaction, other electron carrier (oxidation/reduction) 


18 


5.77E-03 



B. napus genes expressed during infection at 7 dpi were assigned to functional categories and overall pathway regulation was examined. The top 20 categories 
determined by MapMan analysis are listed. Abbreviations: PS (photosynthesis), LHC (light harvesting complex), CHO (carbohydrate), PR (pathogen response), PME 
(pectin methyl esterase), ATP (adenosine triphosphate). AAA type (ATPases Associated with diverse cellular Activities). P values indicate the likelihood of the observed 
pathway regulation being due to chance. 
doi:1 0.1 371 /journal.pone.01 03098.t009 



(Lbjl54_P004138) and Lmb (Lema_P082270.1) (Table 2, Table 
S10, Table Sll in file S2). Its role in disease is unknown. The 
MEROPS SOI A "type" peptidase is bovine chymotrypsin but 
fungal homologs are usually described as "trypsin-like" because 



they are similar to both bovine trypsin, and trypsin from 
Streptococcus bacteria [47] . The in vitro growth condition resulted 
in a more distinct expression profile of peptidases, producing low 
correlation co-efficient values (0.35 and 0.38) between the in vitro 
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and 7 dpi in planta conditions, for Lbc and Lmb, respectively. 
Four Lmb peptidases had moderate levels of expression at 7 dpi 
and were down-regulated during growth in vitro, in a similar 
manner to the trypsin homologs. Two genes, Lema_P058000.1 
and Lema_P03 1600.1, were most similar to C56 cysteine 
peptidases. The other two genes, Lema_P044030. 1 and Le- 
ma_P044810.1, were most similar to the S33 serine peptidases, 
which typically release an N-terminal proline from a peptide 
substrate. These peptidases may target hydroxyproline-rich 
glycoproteins such as extensins, which are components of the 
plant cell wall [48]. Average CAZy class expression for the 
peptidases showed Lbc expressed its peptidase inhibitors more 
highly than Lmb did across all treatments (Table Sll in file S2). 
Glutamyl peptidases were expressed at very low levels in both 
fungi, and threonine peptidases had the highest average expression 
of all peptidase classes. 

Expression of Brassica genes 

The infection stage in a plant-pathogen interaction is often 
reflected in expression of key genes implicated in host defence 
signalling pathways. The expression levels for eight defence 
reporter genes of B. napus were determined during invasion by 
Lmb and Lbc (Figure 5). In general these genes were more highly 
expressed at 7 dpi during infection by Lbc than by Lmb. During 
Lbc infection, reporter genes for ethylene signalling (AGS2, CHI, 
and HEL) were expressed 50- to 138-fold higher than in mock- 
infected controls at 7 dpi, while Lmb infection resulted in only 2— 
19 fold induction of these genes. L. maculans 'brassicae' triggered 
a higher expression level of a key salicylic acid reporters, ICS1 and 
WRKY70, but another salicylic acid reporter gene, PR-1, was not 
as up regulated compared to Lbc. The salicylic acid reporter gene, 
PR-1, but not ICS1, were similarly highly induced at 7 dpi in Lbc. 
The jasmonic acid reporter gene PDF1.2 was induced highly by 
Lbc infection, but not by Lmb. The upregulation of these genes 
reflected that jasmonic acid and salicylic acid defence pathways 
were induced by Lbc, in a pattern consistent with the timing of 
necrosis. 

Sasek et al (2012) showed that interactions of B. napus and 
Lmb, involving the recognition of AvrLml (or the mutant allele, 
avrLml) and corresponding resistance gene, Rlml, salicylic acid 
biosynthesis and transcription of S A- associated genes (ICS1, 
WRKY70 and PR-1) increased as early as 3 dpi. Expression of 
HEL and CHI, genes involved in ethylene signalling, increased at 
7 dpi [49,50]. Although these genes were upregulated during the 
susceptible response compared to the uninoculated controls, they 
were much more highly expressed during a resistance response. 

The increased level of expression of B. napus genes involved in 
ethylene, jasmonic acid and salicylic acid signalling during Lbc 
infections prompted further examination of biochemical pathways 
during infection. A ratio of expression values of B. napus genes at 
7 days after infection with Lmb compared to Lbc was calculated 
and analysed by MapMan software to identify pathways co- 
ordinately responding to early infection of the cotyledon (Figure 6, 
Table 9). The major difference in host response to infection by the 
two fungi was in genes with photosynthesis-associated activities. L. 
biglobosa 'canadensis' infection resulted in massive down-regula- 
tion of genes involved in photosynthesis (e.g. PSI, PSII), electron 
transport, photorespiration and chlorophyll (tetrapyrrole) biosyn- 
thesis compared to that in Lmb. Similarly, expression of sucrose 
and starch biosynthesis genes was reduced, possibly as a flow-on 
from lack of photosynthesis. Levels of sucrose degradation genes 
were higher during Lbc than during Lmb infection, perhaps due to 
decreased amounts of photosynthate. Biosynthetic genes for 
raffinose, a monosaccharide osmoprotectant, were induced by 



Lbc, perhaps due to water stress. In contrast, genes involved in 
starch metabolism (both synthesis and degradation) were tran- 
scribed at high levels during Lmb infection at 7 dpi. The cell wall 
remodelling genes of B. napus that modify fi-glucans, mannans 
and pectin were more highly expressed during infection by Lbc 
than by Lmb [51]. Extensive up regulation of host cell wall 
remodelling genes occurred at 7 dpi as the necrotic lesion was 
formed by Lbc, while Lmb infection had much less impact on 
transcription of these genes. A cohort of genes associated with 
secondary metabolism in B. napus was also plotted using Mapman 
on the same dataset (Figure S3 in file SI). Isoflavone reductase 
genes associated with isoflavonoid biosynthesis, were more highly 
expressed during Lbc infection, while genes associated with 
carotenoid metabolism (phytoene dehydrogenase, zeta-carotene 
desaturase, lycopene cyclases and violaxanthin de-epoxidase) were 
expressed more highly during Lmb infection. These observations 
are in agreement with the increased necrosis observed during Lbc 
infection at the early stages of infection. As well as increased 
expression of secondary metabolism genes during Lbc infection, 
peroxidase, nitrilase, and cytochrome P450 genes were consistent- 
ly upregulated by B. napus at 7 dpi (Figure S3 in file SI). Twenty 
four general peroxidases were upregulated suggesting a strong 
oxidative burst was deployed during infection by Lbc. Forty two 
cytochrome 450 genes, 15 oxidase genes and 34 nitrilase genes 
were upregulated by Lbc infection, which also suggest a strong 
activation of secondary metabolism. 

Summary 

The hemi-biotroph L. maculans 'brassicae' avoids triggering 
host defence during early infection. This is reflected by the finding 
that at seven days post inoculation, L. maculans 'brassicae' 
expresses a large number of genes with no known domains, many 
of them being small secreted proteins. One class of small-secreted 
protein-encoding genes that is highly expressed at this time is 
CBM50 (LysM) genes, which suppress chitin-triggered PAMP 
immunity and evade detection of the fungus by the plant. Also 
avirulence genes are highly upregulated; at seven days post- 
inoculation, two avirulence genes are amongst the top 20 most 
highly upregulated L. maculans 'brassicae' genes in planta. This 
pattern is consistent with the relatively asymptomatic growth 
phase of this fungus at seven days post-inoculation. In contrast, L. 
biglobosa 'canadensis' expresses a high number of cell wall 
degrading CAZy genes during the first seven days of infection, 
consistent with extensive necrosis and a high degree of activation 
of host defence signalling pathways. 

Supporting Information 

File SI Supporting figures. Figure SI, Classifications of 
CBM50 (LysM) domains in L. maculans 'brassicae', L. biglobosa 
'canadensis', Cladosporium fulvum and Zymoseptoria tritici. LysM 
domains from Lmb, Lbc and Z. tritici (formerly Mycosphaerella 
graminicola) with high sequence similarity to ECP6 of C. fulvum 
were aligned using ClustalW. (A) Domains are numbered by 
proximity to N-terminus (#1 is closest). Residues are coloured by 
similarity (black: 100%, dark grey: 80-100%, light grey: 60-80%, 
white: less than 60%). B) Phylogram based on the amino acid 
alignment of LysM domains from panel A. Branch numbers show 
% bootstrap support and scale bar shows amino acid substitutions 
per site. LysM domains assigned to three Positions A, B and C 
based on sequence similarity to the C. fulvum ECP6 sequence. 
LysM domain organisation in ECP6-like predicted proteins in Lm, 
Lb, Z. tritici and C. fulvum are shown in panel C. Mg LysM genes 
are from Z. tritici. Figure S2, Expression of three LysM- 
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containing genes of L. maculans 'brassicae' grown in vitro and in 
planta. Quantitative RT-PCR analysis was performed on RNA 
from Lmb isolate IBCN18 grown in 10% Campbells V8 juice (in 
vitro), and after infection of cotyledons offi. napus cv. Westar at 3, 
7 and 14 days post inoculation (dpi). Expression levels of 
Lm2LysM (Lema_P102640.1), Lm4LysM (Lema_P025400.1), 
Lm5LysM (Lema_P070100.1), were normalised to those of 
gamma-actin (Lema_P099940. 1). Error bars represent one 
standard error of the mean (n = 2-3 biological replicates). Asterisks 
indicate values significantly different from in vitro levels (p<0.05). 
Figure S3, Response of B. napus secondary metabolism and 
large enzyme families to infection by L. maculans 'brassicae' or L. 
biglobosa 'canadensis. 'RNA-seq gene expression values for a B. 
napus unigene set [27] were used to calculate a ratio of expression 
values (log2) for B. napus genes seven days after infection by Lmb 
or Lbc. Ratios were plotted on secondary metabolism gene groups 
using MapMan software on maps for 'Secondary metabolism', and 
'Large enzyme families' maps [39]. A yellow square indicates a B. 
napus gene that is expressed more highly during Lmb infection, 
while a blue square indicates a B. napus gene with higher 
expression during Lbc infection. An expression ratio close to zero 
is shown with a white square and indicates equivalent expression 
during infection by either pathogen. Only genes with expression 
values greater than 10 FPKM were included. MVA is mevalonic 
acid. 
(PDF) 

File S2 Supporting tables. Table SI, Total number of 
RNA-seq reads aligned to reference genomes. Table S2, 
Oligonucleotide Primers. Table S3, Top 100 most highly 
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